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Abstract 

We discuss two cosmological issues in a generic gauge-mediated supersymmetry 
(SUSY)-breaking model, namely the Universe's baryon asymmetry and the gravitino 
^ ' dark-matter density. We show that both problems can be simultaneously solved if 

' there exist extra matter multiplets of a SUSY-invariant mass of the order of the "//- 

term", as suggested in several realistic SUSY grand-unified theories. We propose 
an attractive scenario in which the observed baryon asymmetry is produced in a 
way totally independent of the reheating temperature of inflation without causing 
any cosmological gravitino problem. Furthermore, in a relatively wide parameter 
space, we can also explain the present mass density of cold dark matter by the 
thermal relics of the gravitinos without an adjustment of the reheating temperature 
of inflation. We point out that there is an interesting relation between the baryon 
asymmetry and the dark-matter density. 



On leave from University of Tokyo. 



1 Introduction 



The mediation of supersymmetry (SUSY)-breaking effects to the standard-model (SM) 
sector is a fundamental ingredient in the SUSY extension of the SM. Among various 
models proposed so far, the gauge mediation model with dynamical SUSY breaking 
(GMSB) [in is the most attractive. It may not only solve the flavour-changing neutral- 
current (FCNC) problem, but also provide a natural explanation of the large hierarchy 
between the electroweak scale and the Planck scale because of the dynamical nature of 
SUSY breaking. Extensive analyses of the GMSB models have been performed on the 
physics at the TeV scale. 

The remaining important task is to construct a successful cosmology in the GMSB 
models. The biggest difficulty comes from a stringent upper bound on the reheating tem- 
perature of inflation Tr in order to avoid the overproduction of gravitinos 0. The upper 
bound on Tr is about 10^ GeV when m3/2 — 10 MeV for instance, and it even reaches 
10'^ GeV in the case of the light gravitino m-i/2 ^ 100 keV. In addition, we have to 
adjust the reheating temperature to just below this upper bound if we want to explain 
the required mass density of cold dark matter (CDM) by the gravitino LSPs (lightest 
SUSY particles). Furthermore, we also have to explain the required baryon asymmetry 
with just the same reheating temperature. These facts severely constrain models of infla- 
tion and baryogenesis mechanism, and require a doubtful conspiracy between these two 
independent processes. 

Even only explaining the required baryon asymmetry is very difficult. The unique 
mechanism to generate the baryon asymmetry in such low reheating temperatures, Tr ^ 
10^ GeV, has been considered to be the Affieck-Dine (AD) baryogenesis [| However, 
in recent developments, it becomes clear that the coherent oscillation of the fiat direction 
field (AD field) is unstable against spatial perturbations and that this field fragments 
into non-topological solitons called Q-balls after dozens of its oscillations. In GMSB 
models, it is known that the resultant Q-balls are absolutely stable, and this fact has 
motivated some works to relate the observed baryon asymmetry and the mass density of 
dark matter |]^. Now that some detailed lattice simulations ||^ have been performed, we 
can easily calculate the typical size and the number density of the Q-balls generated in 

^In the case of a heavier gravitino m3/2~ 10 MeV, the leptogenesis through inflaton decays |Q or a 
coherently oscillating right-handed sneutrino |^ can explain the observed baryon asymmetry. 
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GMSB models. Unfortunately, these simulations have revealed that the produced Q-balls 
generally overdose the Universe, which means that the usual AD baryogenesis is not a 
phenomenologically viable scenario as the origin of the observed baryon asymmetry in the 
GMSB models. 

A possible solution to the difficulty is to adopt the AD field, which has a SUSY- 
invariant mass term. This is the case for the AD leptogenesis via the LHu fiat direction 
The scenario can explain the baryon asymmetry provided that the "/i-term" exists even 
in the very high energy (above the intermediate) scale, and that the mass of the lightest 
neutrino is extremely small |jlO|- Q However, we need a conspiracy again between the 
model of infiation and the leptogenesis if we want to explain the gravitino dark matter, 
simultaneously. 

In this paper, we propose a solution to both problems. We consider the cosmology 
in GMSB models under the existence of extra matter multiplets of a SUSY-invariant 
mass of the order of the /i-term. The existence of extra matter multiplets, such as those 
transforming 5 + 5 under the SU(5)gut group, is indicated by several realistic SUSY 
GUT models in which we can naturally solve the doublet-triplet splitting problem 



by utilizing some discrete symmetries. 

As suggested from the AD leptogenesis, adopting an extra matter field as a member of 
the AD fields opens up various possibilities to construct successful AD baryogenesis sce- 
narios. In the present work, we will propose the most attractive one, which can generate 
the required baryon asymmetry in a way totally independent of the reheating temper- 
ature of infiation, without any confiicts with the constraint on cosmological gravitino 
abundance. Surprisingly, we can explain simultaneously, in a relatively wide parameter 
space, the required mass density of CDM by the thermal relics of the gravitino LSPs 
without an adjustment of the reheating temperature of infiation. In this region, we derive 
a novel relation between the mass density of the baryon and that of dark matter, which 
is written solely by low-energy parameters. We will also discuss the non-thermal grav- 
itino production via late-time decays of the unstable Q-balls, and its phenomenological 
implications and experimental constraints. 



^This fact allows us to have definite predictions on the rate of neutrinoless double beta {0v(3(}) decay, 
which will be tested in future experiments |ll[| . 
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2 The Model 



In the present scenario, we assume the existence of a pair of extra matter multiplets, 5x 
and 5a;, which have a SUSY-invariant mass term in the superpotential, W D mx5x5x- f\ 
In the following discussion, we assume that the size of is of the order of the /x-term 
~ 1 TeV), as suggested in realistic SUSY GUT models |T^. 

By adopting one of these extra fields as an AD field, this mass term supplies a quadratic 
potential along the relevant flat direction, which forbids the formation of stable Q-balls 
in the GMSB models. This is not the case, however, when the mass term is generated 
by an expectation value of some field with a renormalizable coupling to the extra matter 
multiplets, such as W D 35x5^, where the S field is to obtain an expectation value. If 
this is the case, the S field is driven to the origin when the AD field develops a large 
amplitude, and then the required quadratic potential vanishes. 

In order to avoid the formation of stable Q-balls in the GMSB models, the mass term 
should be present even after the AD field obtains a large amplitude. This condition can 
be easily satisfied, for example, if the mass term is generated by some non-renormalizable 
operator []T§, W D {S'^ /M^)5x5x, where = 2.4 x 10^® GeV is the reduced Planck scale. 
We assume that this is the case in the remainder of this paper. 

Our model of the AD baryogenesis is basically the same as the one proposed in Section 
II B of Ref. |]T^. We adopt the flat direction labelled by a linear combination of the mono- 
mials of the chiral superfields: UDDx, QD^L. Here, Q, U, D and L denote superfields 
(and their scalar components) of left-handed quark doublet, right-handed up- and down- 
type quarks and left-handed lepton doublet, respectively; belongs to the 5^ multiplet 
and has the same gauge quantum numbers as D. As for an explicit parametrization of 
the flat direction, see Ref. |T^|. 

The relevant baryon-number (B)-violating operators are supplied in the Kahler poten- 
tial. (The effects of possible non-renormalizable operators in the superpotential will be 
discussed later.) We take the following operators, for example, which are consistent with 
the i?-symmetry: 

SK = X.^QUWlL + X.^QU^D^L + h.c. , (1) 



•^The following arguments can also be applied to the case where the extra matter fields belong to 
10 + 10 representations of the SU(5)gut- 
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where Aj's are coupling constants. The SUSY-breaking effect due to the non-zero energy 
density of the Universe induces the following term in the scalar potential, which provides 
the AD field with a phase rotational motion: 

5V = 3H^ 6K , (2) 

where H is the Hubble parameter of the expanding Universe. Here, the fields in SK 
represent the scalar components of corresponding superfields. 

The total scalar potential of the AD field (p can be written as follows: 

V=iml-CHH'M' + ^iaH<l^' + h.c.) + ... , (3) 

where the ellipsis denotes the higher-order terms coming from the Kahler potential; an is 
a complex coupling constant and ch = is a real positive coefficient. The definition 
of (j) is the same as that in Ref. |jl5|, which denotes the complex scalar field along the flat 



direction. The second term is induced from the four-point couplings between the inflaton 
/ and the AD field in the Kahler potential {K D I^I(f)^/M^) 0. 



3 Baryon Asymmetry 

The evolution of the AD field and the resultant baryon asymmetry are almost the same 
as those presented in Refs. |]T^ [T^. Here, we briefly review the important points. The 
AD field is assumed to develop a large expectation value during inflation, because of the 
large negative Hubble-induced mass term. After the inflation ends, the Hubble parameter 
gradually decreases. The AD field starts coherent oscillations when m^. eventually exceeds 
the Hubble parameter. At this time H = Hose — ^x, a huge baryon asymmetry is 
generated by the B- violating operators in the scalar potential (0). 

^ The terms induced by the SUSY-breaking effect in the true vacuum have only a neghgible role in 
the present scenario. Note that the existence of SV in Eq. is independent of the quantum charge of 
the inflaton. This is in clear contrast with the model proposed in Ref. which requires the inflaton / 
to be a singlet in order to write relevant interaction terms, SK = (//M,)(/>t0 and 6W = {I /M,,)WayV°' 
to solve the Q-ball problem in the GMSB models. Such interactions generally lead to the high reheating 
temperature of inflation, which results in the overproduction of the gravitinos. In addition, there is 
no definite description of the available flat direction in the MSSM and the prediction of the resultant 
baryon asymmetry. It is not clear to us whether their scenario can actually produce the observed baryon 
asymmetry without a conflict with the cosmological gravitino problem. 
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One can easily get the ratio of baryon to 0-number density as 





2 

SeS , (4) 



where Ses = sin(arg(a//) +4arg(0o))- This is fixed at H = iJosc and remains constant until 
the fields decay into the gravitinos and SM particles. Here, 0o denotes the expectation 
value of the AD field just after infiation. We treat it as a free parameter for a while, and 
assume that |0|o~ M^. Actually, there exists a natural way to stop the AD field below 
the Planck scale. This can be done by gauging the U(1)b-l symmetry. Although the 
B- violating operators in Eq. (|l]) are invariants of the U{1) b-l symmetry, the relevant fiat 
direction field carries non-zero B — L charges and can be lifted by the D-term potential. 
In this case, the initial amplitude of the field is given by the B — L breaking scale, 
|0|o ~ Vb-L 



In the present scenario, the quadratic potential supplied by the SUSY-invariant mass 
term of the extra matter field forbids the formation of stable Q-balls. However, it 
cannot forbid the formation of unstable ones. This is because radiative corrections of 
SU(3)c and U(l)y gauge interactions give rise to a negative contribution to the mass 
term m^|0p. At the one- loop level, the potential is given by 

^=-^^(i+^iog(|^))i0r, (5) 

where 

1/16 4 \ 

K ^ —Us + — tti . (6) 

47r V 3 15 V ^ ' 

Here, Mq denotes a renormalization scale at which m^. is defined. In this expression, 

we neglect the corrections from the possible Yukawa interactions, which should be small 

enough to avoid large FCNC interactions. In the following discussion, we use K = —0.03 

for a representative value. These radiative corrections make the potential a little bit 

fiatter than the quadratic one, which causes spatial instabilities and the formation of 

unstable Q-balls inevitable. A detailed discussion of the properties of the unstable Q-ball 



can be found in Ref. 14 



^It depends on the matter contents of the model and the scale of \<j)\. However, as we will see, the 
resultant baryon asymmetry and the mass density of dark matter depend only weekly on the precise value 
of the if-factor. 
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For the purpose of calculating the resultant baryon asymmetry, we need the decay 
temperature of the Q-ball. In the present model, the resultant Q-ball charge, Q, is 
estimated as 

where /? ~ 6 x 10^'^ and tc — 0.01 are determined by the detailed lattice simula- 
tions [^, which are independent of |0|o. The size of the charge that evaporates from 
the surface of a single Q-ball through interactions with thermal backgrounds is about 
/S.Q ~ 10^® [|19|. Consequently, the produced Q-balls survive thermal evaporation as 
long as |0|o ^ 10^^ GeV. The remaining charges are emitted through its decay into light 



fermions. The decay rate was calculated in Ref. as 



Here, A = 4:7tRq is the surface area of the Q-ball, where Rq ~ V2/{mxy\K\) is the Q- 
ball radius; u ~ rrix is the effective mass of the Q-ball per 0-number ||2l|. In the present 
model, the above inequality is almost saturated. The detailed discussion of the decay 
process and the condition for this saturation will be found in the next section. 

From Eq. (|), the Q-ball decay temperature can be derived as follows: 
T 

1/2 

QJ • 

where r/ ~ 1 denotes the ambiguity coming from an inequality of the decay rate in Eq. (|^); 
g*{Td) is the total number of effectively massless degrees of freedom at T = T^. From 
Eq. (|^), we can also write in terms of the initial amplitude of the field: 




(10) 



The low decay temperature indicated in Eq. (|^) makes it much easier for the Q-balls 
to dominate the energy density of the Universe. The condition for the Q-ball dominance 
is given by 

Tn>3Tj^] . (11) 
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In the remainder of this paper, we consider the case where the above condition is satisfied. 
The corresponding region in the (|0|o-Tr) plane with thermal effects taken into account is 
presented in Ref. fl^. Basically, there is a wide parameter space for |0|o> IQis GeV. In 
the case of |0|o ~ Mqut 2 x 10^^ GeV), for example, 10^ GeV< Tr< 10^° GeV is the 
viable range of the reheating temperature, which satisfies Eq. (|TT|) . Here, the upper bound 
on the reheating temperature comes from the condition to avoid the early oscillation of 



the AD field by the thermal effects |15 



In the case of the Q-ball dominance, the resultant baryon asymmetry is given by the 
following simple formula: 

/ \ / \ / \ 2 

Pq fnA fnB] 3 



s s \pq J J 4 \M^J 
where s, pg are the entropy density and the energy density of the Q-ball in the Universe, 
respectively. In terms of the density parameter of the baryon number, it is given by 

where h is the present Hubble parameter in units of 100 km sec~^ Mpc~^, and VIb = Pb/ Pc] 
Pb and pc are the energy density of the baryon and the critical density in the present Uni- 
verse. Here, we have expressed the decay temperature of the Q-ball with the initial 
amplitude of the AD field and by using Eq. (0). The corresponding decay temperature 
of the Q-ball can be estimated as ^ 200 MeV for |0|o ^ 10^^ GeV. Note that the re- 
sultant baryon asymmetry is totally independent of the reheating temperature of infiation 
and the detailed history of the Universe , since the final entropy density is solely pro- 
duced via the late-time decays of Q-balls. The observed baryon asymmetry, flBh^ — 0.02, 
suggests the initial amplitude of the AD field to be 10^^ GeV^ |0|o~ 10^''' GeV, which 
surprisingly coincides with the B — L breaking scale suggested from the see-saw mecha- 
nism ra. R 



4 Gravitino Abundance 

So far, we have concentrated on the calculation of the resultant baryon asymmetry in 

a Q-ball-dominated Universe. Now, let us turn our attention to the calculation of the 

^ The mass of the right-handed neutrino Af^j is given by = Xvb-l, where A is some couphng 
constant. This fact and the atmospheric-neutrino data naturaUy suggest vb-l ^ -^gut- 
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relic gravitino abundance. In the present scenario, there are two stages in the gravitino 
production. The first one is the well known thermal gravitino production, which takes 
place at the reheating epoch of inflation. The second one is the non-thermal gravitino 
production from the late-time decays of the Q-balls. 



4.1 Thermal Production 

First, let us estimate the abundance of thermal relics of gravitinos. It was calculated in 
detail in Refs. and is given by 

^3/2 f 45 C(3) /3\ 1 / 90 fm^' 



< ^^^^^^ UJ ' 4 i^^^^j ' ) 

' (14) 
where 

(Sscattt^rcl) ~ 5.9 -^-^ (15) 

is the thermally averaged gravitino production cross section. Here, mg denotes the mass 



of the gluino. The first quantity in Eq. (|Tj) denotes the gravitino abundance when it is 



thermalized; Tc is the "decoupling" temperature of the gravitino, which is roughly given 
by 

C(3) T^(Sscatt^^rcl) 3 



TT 



2 



H{Tc) 2 



(16) 



If the reheating temperature is higher than the mass of messenger (s) quarks Mm, graviti- 
nos are easily thermalized. In this case, Tc should be replaced by the smaller one between 



Mjn and the solution Tc of Eq. (16) 



Thermally produced gravitinos presented in Eq. (0) are significantly diluted by the 
subsequent entropy production associated with the Q-ball decays. The dilution factor via 
the late-time decays of Q-balls is given by 



A T ■ 

^ ini 

where 



Ti,i = Min 



' 3^f^V. (17) 



Tji , Jm^M, 



' 90 

^9* 



1/4 



(18) 
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The second term in the square bracket of Eq. ( |I8D denotes the case where the AD field 
starts coherent oscillation after the reheating process of infiation. In this case, it is difficult 
for the AD field with |0|o — Mqut to escape strong thermal effects, which modify the 
resultant baryon asymmetry and the properties of the produced Q-balls. Therefore, we 
concentrate on the other case, and assume that 

I / qn \^/^ 

TR<^m^M.A^] ^2xl0i°GeV (19) 

in the remainder of this paper. 

First, let us consider the case where the gravitinos are thermalized. From Eqs. (0) 
and (0), we obtain 

s 9*{Tc)Tr \ \(f)\oJ 

In terms of the density parameter, it becomes 

^ ,2 /O.OSX^^V 300 \ / m3/2 \ /lO* GeV\ / \3/2 /iQie GeV\^ 

(21) 

Here, we have used the relation in Eq. (plQl). The required mass density of CDM, ^3/2/1^ — 
0.1-0.2, can be easily explained by adjusting the reheating temperature of inflation. Note 
that this adjustment does not change the resultant baryon asymmetry, since it is inde- 
pendent of the reheating temperature (see Eq. (|13D). 

A much more interesting situation arises if the gravitinos are not thermalized in the 
early Universe. If this is the case, from Eqs. (|l^) and ([17|) , the resultant gravitino abun- 
dance is given by 

One can see that the resultant gravitino abundance is also independent of the reheating 
temperature of inflation. In terms of the density parameter, this is written as 

=^ °- ^ " (^'" ijg" (^) i^r i^r (i^J ■ 

(23) 

Therefore, in this case, both the baryon asymmetry and the mass density of gravitino dark 
matter are solely determined by electroweak-scale parameters and by the initial amplitude 
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of the AD field. Encouragingly, both quantities can be consistent with observations in a 
reasonable parameter space. [] 

In fact, we can derive a novel relation between the mass density of baryons and that 
of the gravitino dark matter in terms of low-energy parameters: 

n 4 / 250 

~ 1.1 X 10 1]'^ 




which is quite consistent with the observation. We should stress that the above combi- 
nation of the baryon asymmetry and the gravitino dark-matter density is independent of 
the unknown high energy physics parameters, such as |0|o and Tr. There remain only a 
weak dependence on the effective degrees of freedom, (^^(Tr), and a natural assumption 
on the coupling constant and the effective CP-violating phase, laj^l^efr = 0(0.01-0.1). 

We can obtain another interesting relation by taking a ratio of the mass density of 
baryon to that of dark matter: 



Q3/2 VlMeVy \ J \ J \ 250 ^ 

As one can see, if the B — L-breaking scale (and hence |0|o) is around the GUT scale, 
the present scenario predicts a successful relation between baryon asymmetry and dark 
matter. 

In Fig. 0, we show a parameter space where we can explain the gravitino CDM inde- 
pendently of the reheating temperature of inflation. (As we have mentioned before, if the 
reheating temperature is higher than the mass of the messenger (s) quarks, the gravitino 



abundance should be calculated by Eq. ([211) 



Finally, we add a brief comment on the possible isocurvature perturbations. In the 
present scenario, the origin of the CDM is the infiaton decay, but the dominant entropy 
is produced by the late-time decays of Q-balls. Thus, the density fluctuation of the CDM 
contains partially isocurvature perturbations if the dominant- density perturbations are 
supplied by the fluctuation of the AD field during inflation [E3]. In order to avoid any 



''Contributions to the mass density of dark matter from the non-thermal gravitino production should 
be subdominant. See the next subsection. 
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Figure 1: The parameter space where the mass density of gravitino CDM is determined 
independently of the reheating temperature of inflation. Below the blue (thick) line, the 
mass density of gravitino is too small to explain the required mass density, Vt^i2h? < 0.05. 
Above the red (thin) line, the gravitinos are thermalized by scattering interactions. The 
region between the two lines, the required mass density of gravitino dark matter, can 
be explained by Eq. (^). Here, we have assumed that the reheating temperature Tr is 
smaller than the messenger (s) quark masses. 



disagreement with observations, the Hubble parameter during inflation should satisfy 
i^/^lO^^|0|o ~ 10^*^^^^ GeV. However, this is not a strong constraint, and it can 
be easily satisfied, for example, in hybrid and new inflation models ||2^. It will be an 
interesting indication of the present model, if a small deviation from the purely adiabatic 



density perturbations is indeed confirmed in the future experiments, such as MAP [25 
and PLANCK 



* In this case, the dominant-density perturbations are generated by the inflaton, and there remain 
small isocurvature perturbations in the baryonic sector. 
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4.2 Non-Thermal Production 



In this subsection, we discuss the non-thermal gravitino production from the late-time 
decays of Q-balls. We first estimate the resultant non-thermal gravitino abundance, and 
then discuss experimental constraints and phenomenological implications. 

As we have mentioned in Section |^, the AD field stored inside the Q-ball decays into its 
superpartner and a gaugino, or a quark/lepton and a higgsino, through gauge or Yukawa 
interactions, respectively. The following condition is necessary for the AD field to decay 
into light fermions with almost saturated rate (77 ~ 1) in Eq. (H): 

/|0|,>^(^m,), (26) 

where / denotes the corresponding gauge or Yukawa coupling constant, and |0|c is the 
amplitude of the AD field at the centre of the Q-ball. We can easily calculate the typical 
amplitude of the AD field inside the Q-ball from its charge and radius, which turns out 
to be |0|c ~ 10^^ GeV for |0|o — Mgut and clearly satisfies the above condition. 

This is not the case for the decay of the AD field into the gravitino. The relevant 
interaction is given by 

2 2 
mz^ — mi _ 

; + h.c. , (27) 

where x<i> denotes the superpartner of the AD field 0, and ip is the longitudinal component 
of the gravitino (~ the Goldstino). In the region where |0| ^ Mm, the mass difference of 
the AD field and its superpartner is given by 

- ^ <ft^ ' (28) 

where msoft — 1 TeV is the soft SUSY-breaking mass of the AD field below the messenger 
scale Mm- As a result, the branching ratio of the AD field that directly decays into 
the gravitino is suppressed by the following factor with respect the decay rate given in 
Eq. (@) 0]: 

<oftMi/\<P\ ^ 1 ^29) 

Below the messenger scale |0| ^ Mm, the numerator in ( P^ is just given by m'^^f^\(j)\. Thus 
the branching ratio is maximized at |0| ~ Mm- The resultant gravitino abundance directly 
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generated from the Q-ball decay is given by 



ubJ Vsm, 



p 



where rrip is the mass of a nucleon. In order to explain the observed baryon asjTumetry, 
the ratio of to Ub is given by {n^/ns) ~ 10^ from Eqs. and (0). Therefore, from 
the above equation, the mass density of gravitinos directly produced by the decays of the 
AD fields inside the Q-balls is a negligible contribution to the total mass density of dark 
matter, as long as Mm ^ 10^° GeV. 

A superparticle produced by the late-time decays of Q-balls subsequently decays into 
the next-to-lightest supersymmetric particle (NLSP x); which is the bino or the lightest 
stau in the minimal GMSB models. The NLSP, in turn, decays into the gravitino and its 
superpartner with the following decay width: 

r ^ (311 

^ " 487rm2/2M2 ' ^ ^ 

where is the mass of the NLSP. In order not to spoil the success of the Big Bang 
nucleosynthesis (BBN), it must satisfy the following constraint [ |28| : 

^ 5Mev( ^-S{U^\>,M.Y. (32) 

VlOO GeV/ V m3/2 / ~ ^ ' 

This condition is satisfied for a gravitino with 7/13/2 ~ (a few) MeV for a reasonable range 
of mass of x- 

First, let us consider the case where 



H{T,) \ = {^^j^\»^.- (33) 



90 M 

Here, the produced NLSPs have enough time to reduce their abundance via annihilations 
into the following value p9[ before they decay into gravitinos: 



n-x I 45 (o-f ) 



-1 



(34) 



13 



where {cfv)^ denotes the s-wave component of the annihilation cross section of the NLSP. 
As a result, the resultant gravitino abundance is given by 

,-2 ( "^3/2 \ GeV-2\ ^100 MeV^^ 



- 2.8 X 10-^ — , " . (35) 

In the case where H{Td) ^T^, the estimation given in Eq. ( pSj ) is not valid, since 
the NLSPs produced by the Q-ball decays do not have enough time to annihilate before 
decaying into gravitinos. In this case, we have to solve the following coupled Boltzmann 



equations fTj]: 



% + 3Hn^ = N^^Q^T' - (^^)x< - , (36) 



^total ^ 3^^total _ Q (f^j, t<Td) , (37) 

n'^""' = (for t>Td), (38) 

^' = ^(PQ+Px + Prad), (39) 

PQ = e-'m^iQ, - TQt)n'^'^' , (40) 
Px = ^x'^x ' (41) 

Prad + 4ffprad = {(^c^'m^ " ^x"^x) Tq^Q*''' + "^x(^^)x^X + "^X^X^ ' (^2) 

n3/2 + = T^n^ , (43) 

where Hq^'^^ is the total number density of Q-balls, which includes positively- and negatively- 
charged Q-balls, ~ m^, is the mass of the AD field, 113/2 is the number density of 
gravitinos, Qi is the initial charge of the Q-ball, = Qi/Tq is the lifetime of the Q-ball, 
and A^^ is the number of NLSPs produced per baryon number, which is roughly equal 
to e~^. Note that the produced gravitinos behave as radiation for the era relevant to 
the present calculation. In Fig. |^, we show the contour plot of figyg/^^ by numerically 
solving the above Boltzmann equations. Here, we have assumed that Td = 100 MeV, 
= 1 TeV, = 100 GeV, e'^ = A^x = 100- 

Now, we discuss experimental constraints on the non-thermal gravitino production. 
First, let us consider the constraint from the BEN. Even if the Q-balls and the NLSPs 
decay into gravitinos well before the start of the BEN (see Eq. (p2D), the produced grav- 
itinos are still ultra-relativistic during the era, which may reduce too much the energy 
density during the formation of light elements. 
^As for the details, see Ref. [Ell and references therein. 
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Figure 2: The contour plot of the mass density of non-thermal gravitinos fig^j^^ the 
(^3/2^(c"^)x) plane. The solid lines correspond to Ogyg/i^ = 0.15, 0.015 and 0.0015 from 
the bottom up. The dashed lines represent the corresponding mass density of gravitinos 
calculated from Eq. (|35D . 

The momentum of the non-thermal gravitino redshifts with temperature: 



where Tg is the cosmic temperature at the time of the gravitino production. Note that Tg 
depends on the production mechanism and the relative size between H{Td) and F^, the 
decay rate of the NLSP, and that Tg is given by 



Therefore, the contribution to the energy density of the radiation from the non-thermal 

^"in the case of the direct production of gravitinos from the AD fields stored in the Q-balls, should 
be replaced by rux- 




(44) 
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gravitinos at temperature T is estimated to 
P3/2(T) ^ f^) s{T)p{T) 




{n^jy)T' , (46) 

where Q^^h'^ denotes the present mass density of non-thermal gravitino dark matter, 
which is nf/^'^^h^ or Vl^/^h^. 

In order not to affect the expansion of the Universe during the BBN, the gravitino 
contribution to the energy density of radiation should satisfy 

^<6N,, (47) 
Pi/ 

where p,^ = ('7r^/30)(7/4)T^ is the energy density of one neutrino species. Agreement 
with observations of light elements requires SN^ = 0.2-1 This leads to the following 
constraint: 

P3/2 ^ ^ 4 ^ ^0-=^ ( ^ ^ ^^^^^ ^""^^ ^""^^ /^^3/>'^ 



< 0.2-1 . (48) 

As can be seen, this constraint is very weak and almost always satisfied for the interested 
region of the parameter space. 

Although the constraint coming from the BBN is very weak, the non-thermal gravitinos 
should not make up a dominant portion of dark matter. This is because they maintain 
relativistic velocity until just before the matter-radiation equality. Actually, they are too 
hot to constitute acceptable warm dark matter, even though they are a little bit cooler 
than neutrinos. The current velocity of the non-thermal gravitino is given by |^ 

( g*iTQ) V'^ To 

2^3/2 \g*{Tg) J Tg 

./lMeV\/ \ /5 MeV\ 

l-SxlO-M^:^ 77^77^ ' (49) 



m3/2 ; VlOO GeVj \ Tg ) 

where Tq = 2.7 K is the current temperature of the cosmic microwave background (CMB). 
Note that the velocity Vq is independent of m3/2 as long as H{Ta) ^ F^, since the gravitino 
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In the case of a direct production, my should be replaced by rrix- 
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mass dependence is cancelled out by the m3/2 dependence of its decay temperature T3/2. 
This is the case when 7713/2 ~ 100 keV for the typical Q-ball decay temperature Td ~ 
100 MeV. In the case of lighter gravitino m3/2 ^ 100 keV, vo is inversely proportional to 
the gravitino mass, since Tg = Td ^ 100 MeV is independent of the gravitino mass (see 
Eq. (^5D). On the other hand, f ^ 1.5 x 10~^ is required for warm dark matter to be 
consistent with observations of the Lyman-a forest 0, which corresponds to the free- 
streaming scale Rf ^ 0.1 Mpc. In the present scenario, Rf Mpc for m3/2 > 100 keV, 
and it even reaches Rf ~ 100 Mpc for m3/2 ~ 1 keV. 

Therefore, the contribution from the non-thermally produced gravitinos should be a 
subdominant component of dark matter; otherwise it would be in conflict with the ob- 
served large-scale structure of the Universe. Although there has been derived no accurate 
bound on the mass density of such "hotter" warm dark matter, we can obtain a conser- 
vative bound by treating it as hot dark matter. CMB measurements and galaxy cluster 
surveys already constrain the mass density of hot dark matter as ^hoth'^ ~ 0.05 at the 95% 
CL. ||33|. From Eq. (|30|) and Fig. ^ one can see that this constraint is satisfied as long as 
M„ < 10^° GeV and {av)^ > (a few) x lO^^^ GeV"^ In the case of the stau NLSP, the 
s-wave annihilation cross section is of the order of 10"'^ GeV~^ (100 GeV/m^)^, which 
clearly satisfies the condition. Even in the case of the bino NLSP, {av)^>10-^^ GeV"^ 
is naturally obtained for relatively large tan/5 (^ 15). If the non-thermal gravitinos make 
up at least few per cent of the total mass density of dark matter, data from the MAP or 



PLANCK satellites may confirm its existence in the near future [34 



5 Conclusions and Discussion 

In this work, we have proposed a solution to the baryon asymmetry and the dark matter 
problems in the GMSB models. We have found that both problems can be simultaneously 
solved if there exist extra matter multiplets of a SUSY-invariant mass of the order of the 
/i-term, which is suggested by several reahstic SUSY GUT models utilizing some discrete 
symmetries to realize the doublet-triplet sphtting in SUSY GUTs |T2|. The resultant 



baryon asymmetry is totally independent of the reheating temperature of inflation, and 
the required initial amplitude of the AD field is perfectly consistent with the B — L 
breaking scale vb-l ~ ^gut suggested from neutrino-oscillation data. Furthermore, 
the abundance of thermal gravitino dark matter is also independent of the reheating 
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temperature in a relatively wide parameter space. In that region, the abundance of 
gravitino dark matter is solely determined by the electroweak-scale parameters and the 
initial amplitude of the AD field. Encouragingly, the required amplitude of the AD field 
is quite consistent with the value required to explain the observed baryon asymmetry. 

We have also discussed the non-thermal gravitino production via the late-time de- 
cays of Q-balls. The constraints from the BBN and the observation of the large-scale 
structures in the Universe can be easily satisfied unless the messenger scale is too high. 
Mm ~ 10^'' GeV, or the gravitino is relatively heavy, ^ 10 MeV. In the case of the 
bino NLSP, the "hotter" warm dark matter supplied by the non-thermal gravitino produc- 
tion is expected to be larger than a few per cent of the total mass density of dark matter, 
whose existence may be tested by MAP or PLANCK experiment in the near future . 

Finally, we add a comment on the effects of possible non-renormalizable operators 
in the superpotential. The discrete symmetry, which gives us a natural solution to the 
doublet-triplet splitting problem, may allow non-renormalizable operators in the super- 
potential, which lift the relevant flat direction. For example, Z^pi symmetry allow 
non-renormalizable operators proportional to (p^, such as oc [UDDy. In the presence of 
the gauged \J{1)b-l symmetry, the operator is written as 6W = (AX/M^)0^, where A is a 
coupling constant, and X is the Higgs field carrying B—L charge 2 to break the B—L sym- 
metry spontaneously. Q This operator lifts the field at |0| ~ M*(ma;/A(X))^/'^, which 
is of the same order as the B — L-breaking scale. Hence, even if this operator exists, it is 
still possible for the (f) field to have its initial amplitude of the order of the B — L-breaking 
scale. The contribution of this operator to the baryon asymmetry is approximately given 
by (nB/n^) ~ (^3/2/"^a')^efT5 which is subdominant as long as ^3/2^ 1 MeV respect to 
that in Eq. (^) (see also the discussion just below Eq. (PDD). Therefore, the results of 
the present model are basically independent of the existence of such non-renormalizable 
operators in the superpotential. 
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